Meson mass spectrum from relativistic equations in configuration space 
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A method is described for solving relativistic quasi-potential equations in configuration space. The 
Blankenbecler-Sugar-Logunov-Tavkhelidze and an equal-time equation, both relativistic covariant 
two-body equations containing the full Dirac structure of positive and negative energy states, are 
studied in detail. These equations are solved for a system of two constituent quarks interacting 
through a potential consisting of a one-gluon exchange part with running coupling constant plus a 
linear confining potential which is mostly scalar and partly vector, and the spectrum of all light and 
heavy mesons is calculated. 



I. INTRODUCTION 



It is rather surprising how well the non-relativistic constituent quark model does describe the masses of the known 
mesons and baryons However, non-relativistic models with or without relativistic corrections can be criticized 
f^f^ ' easily of having an inaccurate treatment of the quark dynamics when high momentum processes or light quarks are 
, considered. It is therefore interesting to study relativistic covariant extensions of the constituent quark model, which 
■ are still likely to reproduce the hadron spectroscopy with similar agreement as the older models, but have a more 
' reliable dynamical framework at extreme momenta. In a previous paper |^ we studied the pi and rho meson in such 
' ' ■ an extension, which was based on a covariant equal-time approximation to the Bethe-Salpeter equation. It was shown 
I that the known high momentum behavior of the electromagnetic form factors could be reproduced within this model. 



. in tnis present worK we investigate two amerent relativistic covariant approximations to tne i^etne-aaipeter equation, 
' namely the equal-time approximation j^] and the Blankenbecler-Sugar-Logunov-Tavkhelidze approximation and 
calculate the meson mass spectrum. The ingredients of the calculations are constituent quarks interacting through a 
^ ' instantaneous flavor-independent phenomenological quark-antiquark potential. The potential consists of a Coulomb- 
^ ,, like one-gluon-exchange part with running coupling constant and of a linear plus constant confining part which is 
^ ' a mixture of scalar and vector character. For the vector part of the potential two different gauges, Feynman and 
Coulomb, are examined. 

The Bethe-Salpeter equation (BSE) is the relativistic covariant generalization of the Lippmann-Schwinger equation 
for a two-quark bound state and can be written as 



S-\P,pMP,p) = - f-^V{P,p-p')i,{P,p'\ (1.1) 

where P and p are the total and relative momenta. Two new features are introduced as compared to the non-relativistic 
equation. The first one is that the bound state wave function ip(P,p) and the potential V{P,p — p') now also depend 
on the fourth component of the relative momentum. The second is that for each quark the number of degrees of 
freedom is doubled due to the introduction of the negative energy states. These extensions make the numerical work 
necessary for solving the BSE highly non trivial. In order to reduce the complexity one may remove all or some of 
the negative energy states, or one may eliminate the dependence on the fourth component of the relative momentum. 
The Lorentz invariance and the structure of the negative energy states can be kept if this component is eliminated. 

There have been various suggestions how to approximate the relative energy dependence of the BSE most of 
which are based on a prescription for modifying the two body propagator in such a way that the relative energy 
is forced to some fixed value. This modification of the propagator is equivalent to a modification of the potential 
and the resulting wave equations are commonly called quasi-potential equations. Here we study two of them. The 
first one is the Blankenbecler-Sugar-Logunov-Tavkhelidze (BSLT) approximation, originally proposed by Logunov 
and Tavkhelidze, and Blankenbecler and Sugar Q, and extended to unequal masses by Cooper and Jennings Q. The 
second is a two-body Dirac equation due to Wallace and Mandelzweig which may be viewed as an extension of 
the equal-time approximation of Salpeter |^ . These quasi-potential equations are normally formulated in momentum 
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space because many expressions in the equations depend on functions like y'pf + rrvf. However, the confining potential 
in momentum space is highly singular at zero relative momentum q. From V{q — 0) — J dx'^ Vi^) and the monotonical 
increasement of the confining potential V{x) one easily finds that V(q ) diverges stronger than q^^ for small momenta. 
Special care must be given when solving such equations I^j9|, |l0| , |ll| . This work uses a different approach. Using 
well-chosen representations of the quasi-potential equations we found that they can be Fourier transformed easily to 
configuration space, thus encompassing all problems associated with the confining potential. 

This paper is organized as follows. In the next section we discuss the derivation of the BSLT and FT equations 
and how they can be Fourier transformed from momentum space to configuration space. The partial wave projection 
of the equations is given in appendix A. In section III we describe the interaction and its implications for the short 
distance behavior of the corresponding solutions. The presence of the Coulomb-like potential makes the wave function 
singular at zero relative distance. This sets an upper bound for fixed coupling constants. The short distance behavior 
of the = 0~ states (e.g. the pion) is analyzed in detail in appendix B. In section IV we describe the method of 
solution of the resulting integro-differential equations. Also we study the parameter dependence and global behavior 
of relativistic effects in the meson mass spectrum. Section V discusses the Regge-trajectories which follow from the 
mass spectra of the BSLT and FT equations. The slopes of these trajectories are found to be dependent on the degree 
of admixture of vector character in the confining potential. If it becomes too large the potential becomes repulsive at 
long distances in some channels making the solutions of the wave equation unbound. This gives an upper bound for 
the degree of the vector character in the confining potential. In section VI the full meson spectrum is discussed as 
obtained from the various quasi-potential equations and in the final section some concluding remarks are made. 



II. THE QUASI-POTENTIAL EQUATIONS 

For the full BSE it is of no importance how the total and relative momenta are defined. However, the prescriptions 
for the quasi- pot ential approximation depend on the definition of the relative momentum. Let us follow Wightman 
and Carding |12| and define the total and relative momenta P and p for a two particle system by 

P = P1+P2, Pi = a{s)P+p, 

(2.1) 

P = P{s)pi- a{s)p2, P2 = P{s)P-p, 

and 

ais) = 1±^±:^, ;,(,) = i^Z!i±Z!i, (2.2) 
Zs 2s 

where s = —P^ = -Pq ^ -f ^'^^ Pi ^^^'^ P2 denote the momenta of the two particles. In the center of mass system, 
where P = (M, o), they read 

Pi = { ^-+po,p) = {E,+pf,,p), (2.3) 

P2 = { ^~po,-p)^iE2-po,-p). (2.4) 

This definition of p has the advantage that in the limit of one infinitely heavy particle, say 7712 — > 00, that Ei — > M—m2. 
Now M — 1712 is precisely the total energy associated with the light particle, so this implies that the relative energy 
Po goes to zero. It turns out that this last property is essential for the quasi-potential approximations to have the 
correct one-body limit, i.e., that they reduce the one-body Dirac equation in the limit of one infinitely heavy particle. 

A. Blankenbecler-Sugar-Logunov-Tavkhelidze equation 

The propagator of the Bethe-Salpeter equation is 

S{pi,P2) [A - "^l)"^(^2 - m2)~^ 

= (A + mi)iA + 5^™'"'^(pi,P2). (2.5) 
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The pole-structure of this propagator is contained in the scalar part of the propagator S'*™'"'"(pi,p2)- In the ap- 
proximation first proposed by Logunov and Tavkhelidze, and Blankenbecler and Sugar it is assumed that the 
pole-structure of the scalar part may be approximated by means of a dispersion relation according to 

^.ca^a,.^^|ca^ar^2^^ / ds' ^ ^ '''^ 5+{[a{s')P' +pY + ml) 6+(\l3{s')P' - + ml) , (2.6) 

where P' — \J s' / sP, and f{s',s) is any function that satisfies f{s,s) = 1. Performing the integral gives in the cm 
system 



^BSLT = -f ^^ 5 (>[P0), (2.7j 

[Ll>i + UJ2) — S IlO\L02 



with uji = -^/p^ + m| . Taking for f{s',s) the simple form proposed by Cooper and Jennings which yields the 
proper one-body limit in the case that one of the masses goes to infinity, one finds 



1 1 2s 



SiPo), (2.8) 



and 

Sbslt{P,p) = (#1 + mi) (#2 + m2) SfslrC^^P)- (2-9) 

The pi and p2 follow from pi and p2 by putting the relative energy to zero, thus pi = (i?i,p) and p2 = {E2,—p). 
Note that this propagator may be cast in the simple forms 

SBSLT = --r ■ -, ^^^'^^ . {fi + TOi)(fe + W2) (2.10) 

[S - [OJI -f CJ2) J [S " [lOi - CJ2) J 2{ijJi + UJ2 ) 



5{vo) ^i+mi 



2(^1 -t- UJ2) #2 - TO2 
^(Po) #2+^2 



(2.11) 
(2.12) 



2{uji + UJ2) f \ - mi 

It is instructive to project the BSLT propagator upon positive and negative energy states. Using 

+ m, = {uJi+poi - is) A+(p,) - (wj - poi - is) KiPi)^ (2.13) 



with 



Ar(p.)- ^-^'^''-^-+;^-^+"'"°' , (2.14) 
and introducing A''^''^ = A^*^ (p)A2^ (— p) gives 

Sbslt = , ^ [(c^i + £^0(^2 + S2)A++ - (c^i + £;i)(c^2 - £^2)A+- 

2(wi -I- W2) (j- 

-(t^i - £;i)(t^2 + -B2)A-+ + {uji - Ei){lu2 - E2)A--] , (2.15) 

where G — uj^ — Ef = cj| — E2- This form clearly shows that the BSLT propagator is time-reversal invariant. It 
satisfies the condition that it does not change form if one interchanges simultaneously the energy labels pi — s- — 
and the particle energies Ei —f —Ei. If the potential is also time-reversal invariant then the BSLT equation has the 
property that if a solution exists with mass M then there will also be a solution at — M, which is identical to the 
solution at M but with positive energy states changed to negative energy states and vice versa. These solutions at 
"M correspond of course to the antiparticles of the solutions at M . 

Since we want to solve the BSLT equation Sggj^rpip = —Vip in configuration space we have to Fourier transform it. 
In order to bring this equation in a form which may be easily transformed we multiply it by 
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D = 



- p2+mi+m2 
2{uji + UJ2) 



which, using eq. (2.11) and (2.12), gives 
(^1+ p2 - nil - m2)i^{P,p) = 



2{lui + UJ2) 



1+ ^2 + mi+ m.2) J -^-^ V{p' - p)iP{P,p'). 



(2.16) 



(2.17) 



In this form it is most clearly visible that in the limit that one of the particles masses goes to infinity, say m2 — > 00, 
the BSLT equation reduces to the Dirac equation 



(#1 - "^i)V'(Pi) 



dp[ 

(27r)3 



m-Pi)V'(p'i). 



(2.18) 



The Fourier transform of the BSLT equation now reads 



nil — 1712 



dx 



, Zbslt{x' - x) 
2(7711 + "12) 



^7(1) . V + £^17^^' - ^(2) . V + £^270 + mi + 7712 V{x')il^{P, x') 



(2.19) 



Due to the relativistic phase space factor in the BSLT propagator, a non-locality occurs in the relativistic equation, 
which is contained in the function 



Zbslt{R) 
This integral can be found by rewriting 



dp Till + ?7l2 Jp . Ji 



(27r)3 uji + UJ2 



Till + 7712 



1 



UJl +UJ2 Till - 7712 

and using the Fourier transform of mi/uji 



UJl UJ2 



d-P Jp ■ R 



Ki{m,R). 



(27r)3 uJi 27r2i? 

Ki is the modified Bessel function of the second kind of order one. The non-locality becomes 

1 1 



(2.20) 



(2.21) 



(2.22) 



^BSLT 



(R) 



™i \ Ko{miR) 



2'n'^R'^ mi — 7712 
2 



Till R 



Ki{miR) \+mil Ko{m2R) 



77l2i? 



Ki{m2R) 



(2.23) 



which is invariant under the interchange mi <-!■ 7712. Both Zbslt{R) and Zi{R) behave like \/R^ for short distances. 
If the two quarks have equal masses then the non-localities become identical, Zbslt{R) = Zi{R) = Z2{R). In order 
to be solved, the BSLT equation is projected on eight basis states with definite parity and total angular momentum 
J. These basis states , th e partial wave analysis and the resulting e ight coupled integro-differential equations are 
discussed in appendix Al, As example, the reader may look at Eq. (B2) to see the explicit BSLT equation for the 
pion. 



B. Equal-time approximation 

In the equal-time (ET) approximation the assumption is made that the qq-potential V does not depend on the 
relative energy in the cm system. This simplifies the BSE considerably. One may then introduce the equal-time wave 
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function = l/(27ri) fdpoip{pQ,p), or il'ix) 

an equation for the ET wave function: 



: i/j{xo = 0, x), and integrate the BSE on both sides over po to obtain 



dpo 
27ri 



S{P1,P2) 



dp 
(2^ 



v{p-p'mp')- 



(2.24) 



For the sake of clarity the derivations in this section are formulated in the cm frame and therefore the ET approach 
does not show a manifestly relativistic covariant form. Yet it can be formulated in fully manifestly relativistic covariant 
way, see Ref. B. 



The ET equation as its stands in Eq. (2.24) does not reduce to the Dirac equation if one of the particles masses 
becomes infinitely heavy and so it does not have the correct one-body limit. However, as has been pointed out by 
Wallace and Mandelzweig |^ , the ET approximation can be improved considerably by adding a propagator Serosa to 
the integral between square brackets, which represents the propagation of the two particles within the crossed box 
diagram. With this extra term the ET propagator exhibits the one-body limit, and furthermore the contributions of 
the crossed diagrams are approximately taken into account. The propagator is now taken as 



Set = J ^ [Sipi,P2) + ScrossiPl,P2)] 



where in the cm system 



ScrossiPl,P2) = S{pi 



cross 
P2 



cross 
P2 



{E2 +Po,-p)- 



(2.25) 



(2.26) 



The only difference between p2 and p^''"'''^ is the sign in front of the relative energy. The origin of the different signs 
for Po can easily be understood from the different directions of flow of relative energy in the uncrossed and crossed 
diagrams. A detailed justification for this crossed box contribution b ased o n the eikonal approximation has been given 
in ref. The integral over po can be evaluated easily by using Eq. (2.13) and expressing the propagators in positive 
and negative energy projections. This gives 



A^ 



A- 



LUl -t- i02 — El — E2 iO\ -f- i02 -|- E\ — E2 



A^ 



A- 



uji + UJ2 — El + E2 uoi + L02 + El + E2' 



(2.27) 



The terms containing A++ and A are well-known and stem from the po-integration over S. They were first used by 
Salpeter 0|. The other two terms result from the po-integration over Scross- This form shows that the ET propagator 
is time-reversal invariant, similarly as the BSLT propagator. The most important difference between the BSLT and 
ET approximation is the degree to which the propagation of negative energy states is suppressed as compared to the 
propagation of the positive energy states. For particles of equal mass {lji = lu2 = uj, and Ei = E2 = E), the ratio of 
the propagation for positive and negative states is for the ET case 



whereas Eq. (2.15) gives for the BSLT case 



S 



ET 



■'ET 



E 



E' 



(2.28) 



S 



BSLT 



^BSLT 



E 



E 



(2.29) 



Therefore we can expect a larger admixture of negative states in the ET solutions than in the BSLT solutions. Note 
that this ratio is always positive for the BSLT propagator whereas it can be negative for the ET propagator. For the 
full BSE it is always positive at po = 0. 



Writing out the projection operators A''^''^ gives 



c-l 



1 

UJl 



mi 



7*^^ ■p-lo^^E2+m2 



1 

U!2 



-l^'^-p-l'i^Ei 



mi 



7(2). p- 



7712 



(2.30) 
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The ET bound state wave equation S'^^V' = —V^ may now be transformed easily to configuration space and becomes 



dx' 



with 



Zi{x - x') Z2{x-x') ■ 
Si{x ) H S2ix ) 



mi 



m2 



Si{x) 
S2{x) 



17 



(1) 



mi 



(2) 



^{x') = -V{x)^{x), 



i^i) . V -7^^^£;i +mi 



-17 



(2) 



m2 



(2.31) 

(2.32) 
(2.33) 



and the non-locaUties Zi have already been given in Eq. ( 2.22 ). The ET equation can be solved after projecting it 
on basis states with definite parity and total angular momentum J. The resulting eight coupled integro-differential 
equations can be found in appendix A 2 . 



III. INTERACTION 



It is commonly accepted that the interaction between the two quarks consists of a short-range part describing the 
one-gluon-exchange (OGE) potential and a infinitely rising long-range part responsible for the confinement of the 
quarks |lJ. We use 



V{x) = -^7(i)7A'(2) + (^^ + c) [(1 - £)l(i)l(2) + £7(1)7^(2) 



(3.1) 



(The color factor 4/3 of the expectation value of the OGE in the meson color wave function has been absorbed in 
the definition of a.) The OGE potential is a pure vector interaction. The confining potential is commonly believed 
to be purely scalar. However, we choose a confining potential which is mainly scalar, but it can have a fraction e 
of vector confinement. Asymptotic freedom requires that for short distances the OGE coupling constant decreases 
logarithmically as a{x) ~ ao/hi{xo/x) where ao = 87r/(33 — 2rti?) and xq — e^'^ / Aqcd |13|- Richardson ||l^ has given 
an elegant prescription for a running coupling constant which for small distances reproduces the correct asymptotic 
freedom and which gives for large distances a linearly rising potential. We do not use this prescription since it can 
not specify the vector and scalar character of the potential. At large distances the confining interaction dominates 
and the exact behavior of a(x) becomes of little importance. If one assumes that the coupling constant grows to some 
saturation value agat then a smooth interpolation between the short and long range is given by |13] 



ai{x) = ao In ^ 



Xo 



exp(— /j,x) -I- exp 



asat 



(3.2) 



Here the typical range of the running coupling regime is controlled by the parameter /i; the maximum range is found 
at fi — 0. A longer range can be obtained by interpolating 



aiiix) = ao 



In 



ao 



asat 



(3.3) 



We use Aqcd — 0.2 GeV and np ~ 3. The dependence on Aqcd and np, which is not large, can be compensated 
for by modifying fi and asat- 

In the BSE the solutions are independent of the choice of gauge for the vector part of the potential. This indepen- 
dence depends essentially on the presence of crossed diagrams in the interaction kernel ||l^. However, in the BSLT 
and ET approximations the relative energy is fixed, so crossed diagrams do not occur and the gauge-independence of 
these quasi-potential equations must be broken. In order to e stim ate the gauge-dependence of our results we study 
two gauges. The first one is the Feynman gauge as used in Eq. ( |3.lD . The other is the Coulomb, transverse or radiation 
gauge which is obtained by replacing 



.,a) ./2) (7W.a.)(7(^)-: 



(3.4) 



in Eq. (p.l|). Appendix |A 3| gives the partial wave projections of these potentials. 
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The presence of the 1/x term in the potential has important consequences for the short distance behavior of the 
wave functions. This can be understood from the foUowing simple picture. Consider the probability of the meson 
decaying through the annihilation of the quark and the antiquark. This process is proportional to the wave function 
at zero relative distance, or equivalently, to the wave function integrated over all relative momenta. Using the wave 
equation the integral over the wave function can be expressed as a J dpip{p) — J dpS{p) J dqV(j> — q)'4'{q), that is, a 
loop integral over the propagator S and potential V plus additional corrections. This is illustrated diagrammatically 
in Fig. ^ Let us focus on the ultra-violet (uv) behavior of the momentum integration over the triangle diagram. It is 
divergent for most dynamical models due to the OGE potential. For example, in the BSE the two-fermion propagator 
and the OGE interaction both fall off as for large p (we neglect for the moment the additional Inp^ behavior of 
the running coupling constant). So for large relative momenta the loop integral takes the uv-divergent form ^ d'^p p~^. 
Similarly, in the one-body Dirac equation and in the two-body quasi-potential equations the propagators behave as 
p^^ and the potential as p~^ what also leads to a uv-divergent loop / d'^p p^^. Divergences also occur in the light-cone 
formalism where the propagator and potential go as p']_^ and the two quark spinors as p\ | [l7|] resulting in Jd^p± p±^ ■ 
This suggests that the solutions of these wave equations are singular at the origin. Indeed, at short distances the 
Dirac wave functions and the Bethe-Salpeter wave functions p^ , pO| behave like "tjjix) '--^ x'^ , — 1 < 7 < 0, where 7 
is a decreasing function of the coupling constant a of the interaction. At some maximum value of a the exponent 7 
becomes smaller than —1 and physical acceptable solutions do no longer exist. These features are also shared by the 
BSLT and FT equations. 

The extent to which these singularities actually appear in the physical wave functions is weakened, if the fixed 
coupling constant is replaced by a running coupling constant such that a{p) ^ In^^ p^ for p — > 00. The triangle loop, 
however, is still uv divergent but the short distance behavior is like ~ \^\- Furthermore, the singular behavior 
can be removed by renormalizing the wave function by means of some cut-off scale. For example, this has been done 
in the light-cone calculations of Ref. |17|. Yet, the mathematical implications are important. Seemingly no attention 
has been given to the asymptotic behavior in most work on relativistic quark-quark dynamics such as in Refs. pl| , p2| . 
In appendix B we analyze in detail the singular behavior of the wave functions of the = 0~ states (e.g. the pion) 
in the BSLT equation with a fixed coupling constant. We have already discussed the short distance behavior for the 
FT equation in Ref. |]. 



IV. PARAMETER DEPENDENCE 



In this section we describe the calculational procedure to solve the quasi-potential equations and we present some 
results on the parameter dependence of the spectrum of J^'~' = states. For convenience we refer to these states as 
the ^5*0 states since their wave functions contain mostly ^So components. Similarly, other J^^ states are also named 
by their main ^^~^^Lj components. Numerically stable solutions are obtained by taking explicitly into account the 
singularity in the bound state wave equation at x = due to the presence of the OGE term. In appendix B the behavior 
of the wave function of the ^Sq state is analyzed in det ail for sma ll x. W e fin d tha t the most singular component of it 
behaves as V' ~ a:"^, -1 < 7 < 0, and 7 given by Eqs. (^), ( ^ ), ( ^TT| ) or (^ ). Let us introduce a function f-y{x) 
which for small arguments behaves as x'^ and for large arguments becomes unity. By substituting ^^{x) — f-y{x)ip{x) 
we get an equation for ip{x) which is regular at x = 0. For the non states, no special precautions have to be taken, 
because the wave function for the coupling constant strengths considered in this work vanishes sufficiently fast in the 
origin. Numerically accurate solutions are found by using the same x'^ factor as for channel. The same applies 
when we have a running coupling constant, since the equations in this case are less singular than those with a fixed 
coupling constant. 

The partial wave projection of the BSLT and FT equations leads to the coupled set of differential-integral equa- 
tions (A9) and ( A20 ) for n (3<n<8) components ipi{x), which are expanded as tpi{x) = f-y{x)^'^^-^CijSj{x). The 
Sj(x) are cubic Hermite spline functions (see e.g. Ref. psf). The ipi{x) are cut off at some maximum value Xmax and 
we impose ipi{xmax) = '^'[{xmax) = 0. Near x = V'i(a^) is forced to be of order x'^ or higher. By evaluating the equa- 
tions ( A9) or (A20) at k fixed points xi, Xfc one obtains a set of fc x n linear equations for the kxn spline coefficients 
Cij, which only admits non-trivial solutions at the bound state energies. It should be noted that due to the multiplica- 
tion of th e BSLT equation by the operator D of Eq. (2.16) a continuum of additional solutions have been introduced in 
Eq. ( 2.17| ) for M > nii +m2, which are not present in the original BSLT equation. The cut-off makes these continuum 
unphysical solutions into a discrete set, thus rendering it possible to isolate and reject them. This is illustrated in 
Figs, ga and b which show two ^Sq solutions of the BSLT equation with m = 0.2 GeV, n — 0.2 GeV'^, e ~ a = c = 0. 
The wave function shown in ^ is the ground state found at M = 1.23 GeV; the solution shown in |^ found at 
M = 1.30 GeV can easily be be identified as an outgoing wave. As a result the latter solution should be rejected. The 
masses found for the physical bound states are insensitive to variations of x^ax , whereas the masses of the continuum 
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solutions do depend on Xmax- Typical values that we used are Xmax = 2 fm and fc = 30 for the heavy quark systems 
up to Xmax = 5 fm and fc = 120 for the pion system. For the latter system a large number of splines is needed to 
obtain good accuracy since its wave function can have a long oscillating tail. As an overall check on the method the 
configuration space program was also used to calculate the spectrum for a non-confining potential containing scalar 
and vector exchanges. The resulting masses were verified by solving the BSLT equation in momentum space [ p4[ ; 
agreement was found within 0.1%. 

Let us now consider some important differences between the non-relativistic and the relativistic equations. In the 
conventional non-relativistic models as discussed for example in Ref. [Q the hyperfine interaction resulting from the 
OGE contribution is singular near the origin and as result it has to be regularized by introducing a phenomenological 
cutoff. On the other hand, our relativistic equations are well defined, at least up to a critical coupling constant of 
the OGE interaction. In this case a natural cutoff essentially occurs due to relativity, where the scale is given by the 
mass of the constituent quarks. The singular behavior of the wave function at a: = induced by the Coulomb-like 
interaction in the relativistic equations considerably modifies the meson mass if the the coupling constant a is close to 
its maximum value. Fig. || shows the masses of a mostly non-relativistic ^5*0 system {m — 5 GeV and k = 0.2 GeV^) 
as a function of the coupling constant. For small a the spectrum agrees well with the Schrodinger result. However, if 
the limit a Umax is taken, the relativistic results differ considerably; dM{a)/da —oo and M{a) Mq, where 
generally Mq is not equal to zero. 

Furthermore, for the non-relativistic Schrodinger-like equations an additional constant c in the potential causes 
solely an overall shift in the meson mass spectrum. This does not happen in the relativistic case. Because the 
quasi-potential equations are time-reversal invariant each solution for a meson at some positive mass is accompanied 
by its anti-meson solution at the same negative mass. If the constant c is added to the potential and we let c grow 
to negative values, the absolute mass of both solutions decrease until both mesons become massless. At this c the 
solutions coincide and for more negative c no bound state can be found except for higher excitations. This can clearly 
be seen from Fig. ^ where the ^5*0 spectrum is plotted as a function of c for the BSLT and ET equations. 

The replacement of the non-relativistic kinetic energy /2m by the relativistic expression y/p^ -\- m? and the 
introduction of the negative energy states greatly reduces the spacings between the ground state solution and the 
successive excitations for small quark masses. Fig. || shows the calculated ^5*0 mass spectrum as a function of the 
quark mass with the confinement taken such that njw? is fixed. This confinement gives for the Schrodinger equation 
constant binding energies. The figure illustrates nicely that for large masses the non-relativistic results are obtained, 
whereas for small masses the level density becomes considerably higher. 

Finally, another aspect which is absent in non-relativistic models are the solutions describing a heavy quark and a 
light quark-hole. Their presence is a consequence of the requirement that the one-body Dirac equation is obtained in 
the limit that one particle becomes very heavy, and they can be interpretated in an analogous way as the negative 
energy states in the single fermion Dirac hole theory. Fig. ^ shows the ^5*0 mass spectrum for fixed mass of the first 
quark and various masses of the second quark. At equal masses the hole-quark states are not present, but as m2 
grows they emerge from the zero total mass axis, and at large two spectra symmetrical around M — ±m2 are 
generated. 



V. LONG DISTANCE BEHAVIOR AND REGGE-TRAJECTORIES 



Experimentally it is found that the masses of the light mesons lie on a linear Regge-trajectory, that is, for large 
angular momenta J the squares of the masses M of the mesons are proportional to their angular momenta 



(5.1) 



where the Regge-slope /3 ~ 1.2 GeV^ ||T^. The value of the slope depends almost exclusively on the confinement 
strength and can be used to fix it. Figs^ and ^ show the Regge behavior as obtained from the numerical solutions 
of the BSLT and ET equations using ruq = 0.25 GcV, a — Q^c— —1.0 GeV and a linear scalar confining potential 
with K = 0.33 GeV^ The fi gures indeed show a linear relation between and J. However, the Regge-slope is 
considerably smaller than predicted by non-relativistic models with energy operator ^/p^ -\- m? , where f3 = 8k, [ p^ ; 
these models give a good description of the cc and bb systems when k ~ 0.18 GeV^ (see e.g. 0). It turns out that 
a small admixture of vector confinement e, as in Eq. (^), greatly affects the Regge-slopes. Figs. ^ and ^ also show 
the Regge behavior for a fraction e = 0.15 of vector confinement in the Feynman gauge and in the Coulomb gauge. 
All trajectories are increased except for the spin-singlet trajectory in the ET model in the Feynman gauge which is 
decreased. Table | summarizes the resulting Regge-slopes. 

A vector admixture in the confining interaction introduces a spin dependence through the presence of the | — h) , |H — ) 
and I ) components. Because of this when too much vector confinement is chosen, e > 0.2, some light mesons are no 
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longer bound, i.e. certain channels may become deconfined. This can be used to set an upper bound on the degree of 
allowed vector admixture. Consider the strength of the confining potential in the ^ Jj channel in the Feynman gauge. 

From Eq. (A22) we see that the projection of this channel is —4, so the linear part of the qq-interaction 



in this channel is 



(1 - e) - 4e 



(5.2) 



Hence if e > 0.2 the potential becomes repulsive for large distances; the interaction tends to — oo for x oo, which 
is clearly not physically admissionable. This suggests that the P = (— mesons of the light quark system, of which 

the wave functions contain the ^ Jj channel and have | ) components of size comparable to the |+ +) components, 

do not have a bound state solution if e > 0.2. This is confirmed by the numerical solutions. In the Coulomb gauge a 
slightly larger value is allowed, namely s — 0.25. In order to avoid these problems we demand that the admixture of 
the vector interaction in the confining potential does not reach or exceed these values. 

It is interesting to note that also in the one-body Dirac equation, which can be used to describe a system of 
a very heavy and a light quark a similar problem arise. Such a framework has recently been used to study the 
spectrum of the D and B mesons p5[ |. A pure scalar confinement interaction was used. The repulsion exercised by 
the vector confinement on the negative components gives also in this case an upper limit on the allowed degree of 
vector admixture in the confining potential. Consider the Dirac equation 



(m — i ^)ip{x) — —nx[{l — e) + e^o\ip{x). 
Projection upon states of definite J gives 



m 
D 



E D, 



J+3/2 

- ,7+1/2 m + E 



1 

1 - 2e 



tp{x), 



(5.3) 



(5.4) 



and Dg as defined in Eq. (A3). Eliminating the lower component of tpi^) iii favor of the upper component of 
tpi{x) — F[x)/x and letting a; ^ oo leads to 



dx"^ 



-|- + 4(-^-3/2)(J+l/2) 
x ox X 



+ (1 - 2e)n^x^ \ F{x) = 



(5.5) 



For e > 0.5 the potential term becomes negative and the solutions for F{x) are no longer bound. 



VI. MESON MASS SPECTRUM 

In Table || the parameters are listed that we used to calculate the meson mass spectrum. We did not include 
the light isoscalar mesons (77, /, h). The masses of these mesons are expected to be modified by the process of qq 
annihilation and creation, making the transitions uu«->dd<-!-ss likely. 

We determine the strength of the confinement from the condition that the experimental Regge-slope /3 ~ 1.2 GeV^ 
is found. If the confinement is taken purely scalar then a confinement is needed which is much stronger than the 
strength known from non-relativistic models to give a good description of the cc and bb systems. Therefore the 
maximum of allowed vector admixture is chosen which gives a maximum Regge-slope at a fixed confinement strength. 
So in the BSLT model we take e = 0.2 and e = 0.25 in the Feynman gauge and in the Coulomb gauge, respectively. 
With these values the experimental Regge-slopes are found at k ~ 0.33 GeV^. Although this strong confinement 
makes it difficult to get a perfect fit for the heavy quark systems we need it in order to get an acceptable light meson 
description. As can be seen from Table | and Figs. and ^, an increase of vector confinement in the Feynman gauge 
in the FT model gives an increase in the Regge-slopes at the spin-triplet states but a decrease at the spin-singlet 
states. It is therefore impossible in the FT model in the Feynman gauge to obtain simultaneously a fair description 
of the heavy mesons and of the Regge-slopes of the light mesons, so for the fit of the meson mass spectrum we only 
used the FT model in the Coulomb gauge. For this last case we took the value e = 0.2 which is slightly less than the 
maximum e = 0.25; above e ~ 0.2 the spatial extension of the ^ Jj components becomes more than '--^lO fm which is 
unphysically large. 

For the OGE part of the potential the choice of the saturation value a^at of the running coupling constant is not 
important since it only plays a role at large distances were the confining part dominates. Fair fits can be found 
with asat = 0.6 as well with a^at = 1-0. We choose the commonly accepted value asat = 0.8. The typical range of 
the running coupling constant regime is much more important. In the Coulomb gauge the maximum fixed coupling 
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constants allow ed in the BSLT and in ET model are respectively amax = 8/(37r) ~ 0.85 and amax = 4/(37r) ~ 0.42 



[see Eqs. ( Bll ) and (B12)]. If the running coupling constant regime is small and the saturation value is approximately 
equal or larger than the maximum allowed fixed coupling constant, asat^ctmax, then at medium small distances the 
wave function exhibits the singular behavior of the fixed coupling constant equations. Especially, the spectrum shows 
much lower energies for the ^5*0 states than the non-relativistic spectrum. Since amax is rather low in the ET model 
this effect occurs rather strongly in the ET spectrum, and one must choose a much larger running couplin g re gime 
for the ET equations than for the BSLT eq uati ons. So we use for the BSLT model the coupling aj of Eq. (|3.2[) and 



for the ET model the coupling an of Eq. (3^). Yet, even with this choice the ^5*0 states are considerably lower in 
the ET fit than in the BSLT fits. 

The masses of the u and d quarks are taken equal and chosen together with the constant c to give a fair description 
of the light non-strange mesons. Next the masses of the s, c, and b quarks are fitted to the '^Si states K* , J/tp and 
T, except for the c mass in the ET model where a correct K* mass would lead to a far too low K mass. 



Table HI presents the resulting mass spectra of the BSLT approach in both the Feynman and in the Coulomb gauge 
and of the ET approach in the Coulomb gauge, together with the known experimental mass spectrum. As can be 
seen from these numbers the BSLT approach in the Feynman gauge gives the best description, but the differences 
between the BSLT results obtained from the Feynman gauge and those from the Coulomb gauge are minor. The ET 
model does not give a very good spectrum. 

Let us discuss the spectra in more detail. The ground states of the dd and ds systems, the n and K mesons, are 
considerably lighter than the BSLT fits predict. The light masses of these mesons are commonly explained within in 
the framework of broken chiral symmetry where these mesons correspond to almost massless Goldstone bosons. Since 
the quasi-potential equations do not incorporate the chiral symmetry they give too heavy tt and K. The correct mass 
of the TT in the ET is a coincidence and due to the too singular behavior of the ^5*0 states in the ET model. 

The description of the ^Pi-states in the ET model is not good. According to the Breit-reduction of the potential, 
the mass of this state should be equal to the center of gravity of the states up to order l/m^. Clearly none of the 
"^Pi-states in the ET model satisfy this condition, in contrast to the BSLT model where this condition is more or less 
satisfied. Furthermore the splittings between the various ^Si and '^Di states, especially in the dd and ds system, are 
much too small in the ET model. We conclude that the ET model does not give the correct fine structure. At this 
point we would like to note that within a constituent quark model with flavor independent potential it is not possible 
to have simultaneously a correct description of the P-states in the dd and ds systems. The only difference between 
these systems is that one of the constituent quarks is a little heavier, so it can be expected that the levels of the ds 
system are raised a little as compared to the levels of the dd system, with a little smaller fine-structure splitting. But 
the observed spacings between the P-states of the dd and ds mesons do not follow this pattern. 

Considerable deviations between the experimental masses and the calculated masses appear in the excited P-states 
of the bb system. This is due to the large confinement strength k that was chosen in order to have proper Regge- 
slopes for the light mesons. We have searched for various simple prescriptions to improve on the shape of the potential 
in such a way that the masses of the excited P-states of the bottonium system show better agreement. The large 
confining potential needed to get the proper Regge-slope for the light mesons and the weaker confinement needed for 
the heavy excited P-states in bottonium clearly suggest that maybe a better overall description can be found if one 
takes a smaller k at short distances than at large distances. We examined this possibility by modifying the potential 
at short distances according to 

{k'x + c' + Ax'^ X < Xo, 
(6.1) 
KX + C X > Xo, 

with the two regions smoothly matching, k' < k, A > and various character (vector or scalar) for Vconf — Vconf- 
We found that it was indeed possible to get the right Regge-slopes for the light mesons and simultaneously the bb 
excitations at the correct levels. But this also increased the fine-structure splittings in the non-strange P-states up to 
values as 0.3 GeV for the difference M{a2) — M{ai) which is known to be only 0.051 GeV. Furthermore, the decrease 
of the splitting M(2^Pi) — M(l^S'i) in bb causes almost a similar reduction of the shift Af(l^Pi) — M{l'^Si) in cc 
which is undesirable. These two differences depend more or less in the same way on the potential since the excited 
levels 2'^Pi of bb and I'^Pi of cc both have approximately the same spatial extension. 

VII. CONCLUDING REMARKS 

We studied two approximations to the Bcthe-Salpeter equation for the two quark system, the Blankenbecler-Sugar- 
Logunov-Tavkhclidze (BSLT) equation and an equal-time (ET) equation. In these quasi-potential approximations 
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the relative energy dependence of the wave function was eliminated by assuming a simplified form for the two-body 
propagator. The full Dirac structure of positive and negative energy states was kept. Both two-body propagators 
reduce to the one-body Dirac propagator if one of the particles is taken infinitely heavy. We applied these equations 
to a system of two constituent quarks interacting through a phenomenological potential which consisted of two parts. 
The first part was a Coulomb-like one-gluon exchange (OGE) part, the second part was a linear confining potential 
which was taken mostly scalar-like and partly vector-like. For the vector part we studied both the Feynman and the 
Coulomb gauge. Since the confining potential is highly singular in momentum space we transformed the quasi-potential 
equations to configuration space. 

It was shown that for a fixed coupling constant for the OGE potential a maximum value exists — depending on the 
model and gauge — above which the ground state solutions no longer exist. Also the admixture of vector character 
in the confining potential has a maximum value (20% and 25% in the Feynman and Coulomb gauge respectively) 
above which some mesons become unbound. The latter could be explained from the repulsion between positive and 
negative energy states. We found linear Regge-trajectories for both models and gauges, but their Regge-slopes were 
much smaller than predicted by models with only positive energy states. 

Using the BSLT and FT equations we calculated the full known meson mass spectrum of all light and heavy mesons. 
As compared to the non-relativistic model predictions [Q given the limited number of parameters used here, the fit 
can be considered satisfactory. We believe that our prediction of the spectrum is of comparable or better quality 
than other relativistic studies based on quasi-potential and Dirac equations |^,^. We found only a small gauge 
dependence in the BSLT spectrum. For the FT equation we only used the Coulomb gauge since it was impossible to 
get satisfying Regge-slopes in the Feynman gauge. The fine-structure of the spectrum, such as the spacings between 
the P-states and between the ^Si and ^Di states, as calculated from the FT equations did not follow the experimental 
spacings nor did it follow predictions of the Breit-reduction of the potential. The BSLT fine-structure showed more 
or less agreement with these. However, from the agreement that was found between the calculated meson mass 
spectrum and the experimental spectrum we conclude that for a calculation able to reproduce all meson masses 
within ^ 0.03 GeV one cannot suffice with a constituent quark model and a phenomenological flavor independent 
potential. Specifically, one would like to incorporate chiral symmetry in order to reproduce the correct tt and K mass 
and some flavor dependence in the confinement strength in order to reduce this strength in the heavy quark systems. 

The wave functions found in these calculations can in principle be used to perform relativistic covariant calculations 
on interactions between mesons and other particles. Especially at high momenta significant deviations from the non- 
relativistic calculations can be expected; one example has been discussed in Ref. M. 
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APPENDIX A: PARTIAL WAVE PROJECTION 



In this appendix we describe the projection of the potential and the BSLT and ET equations on sets of basis 
states with definite parity and total angular momentum. The equations and basis states are formulated for a two 
quark system; they can easily be transformed to a quark-antiquark system by performing a charge- conjugation on 

one of the quarks. The complete meson wave function is written as a combination of the states | JLS'Jz> (E)|pi^''/9i^'> 
^Xcoior'^X flavor- The color and flavor parts of the wave function are well-known and can be found e.g. in Ref. [|l|. For 
the p-spin basis states we do not use the energy eigenstates projected out by the A^^p^ , but instead the four eigenstates 
of 7o^^7o^'', corresponding to the four combinations of upper and lower components of the quark spinors. T he en ergy 
eigenstates have the advantage that the propagator is simply diagonal on this basis as is shown by Eqs. ( 2.15 ) and 
(2.27). However, on this basis the matrix elements of the potential are complicated and many in number. In principle 
it is possible to transform them all to configuration space, but the numerical work involved is huge. On the other 

hand, on the basis of the eigenstates of 7o^^7o^'' the matrix elements of the potential are almost trivial, whereas the 
projection of the propagator on this set gives only a limited number of functions to be transformed. For the BSLT 
propagator there is only one such function namely Zbslt (2.23), whereas for the ET propagator two functions Zi 
and Z2 (2.22) are needed. 
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1. Blankenbecler-Sugar-Logunov-Tavkhelidze equation 



We decompose the kinetic part of the BSLT propagator into pieces acting in the subspaces of the spin and rho-spin: 



(2) 



(2) 



2 Pa ■ 



(Al) 



For a given total spin J there are four basis states for the spin subspace. In the spectroscopic notation "^^^^Lj they 
are {^Jj, ^Jj, ^{J—l)j, ^{J . After some Clebsch-Gordan algebra one finds that on this basis 



r(l) . V = 




(A2) 



where 



ox X 



ci = 



J 



2J+ 1' 



C2 = 



J+ 1 

2J + 1' 



(A3) 



The matrix of cr^^) • V is identical, except for the first row and first column — these correspond to the antisymmetric 
spin singlet state — which change sign. Let us combine the four eigenstates of 7o^''7o^'' into the rho-spin combinations 



l + +>-|--> 



l ++>+|--> 



V2 ' V2 ' ' 

Phase factors i between the eigenstates are used according to 



V2 ■ 



!+-)-!-+) 

V2 



-+) = -lp\ 



■ (2) I 

= -Wi I 



(1) (2) I 
= Pi Pi I 



(A4) 



(A5) 



This ensures that the final equations will be real. Following Kubis and Gammel ||3^] we combine the spin and 
rho-spin bases to form the basis set B of eight states with unnatural parity P = (— )'^^^ 



6 = {Vj, ij^, V}, '(j+im, 

and the basis set B* of eight states with natural parity P — {^Y 

V}, \rj, ^(j-i)}, '(J-m, '(J+i)}, '(J+ir/}. 



(A6) 
(A7) 



These basis states for two quarks can be related to the often used Dirac matrix set for a quark and an anti-quark |3^ ] 
by performing a charge conjugation on one of the spinors. The relation is shown in Table Let us expand the 
expression for f)2 on the set B (A6) of unnatural parities and denote the resulting 8x8 matrix by D{x). This 
gives 

\ 



(A8) 
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■1+1 





C2DJ+2 
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D 
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C2D-J^ 


-1 -ci 


Dj+2 


y ~ 




~ciDj+i 
















A 
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-C2DJ+1 


A 

















C2D-J 






















A 




V 





















A 





+ E2 = 




M and 2A = 


= El 


-E2 = 


{mj - mD/M 


The matrix D* oi p 


L+ P2 



natural parity states can be obtained from D by interchanging E' A. The complete partial wave projected BSLT 
equation is the following set of eight coupled integro-differential equations 



[D{x) - M]ip{x) = 



1 



2(toi -|- TO2) 



dx' Z 



BSLT 



{x,x')[D{x') + M]V{x')ip{x') 



(A9) 
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The matrix M is equal to mi + m2 times the identity, and the diagonal matrix Zbs lt has as i-th diagonal element 
the projection of Zbslt on the orbital angular momentum £i of the i-th basis state, 

zfsLT^^.^') = 27ry d7 F,(7) Zbslt{.^J ^ x'^ - '^xx'-{). (AlO) 

The number of coupled channels becomes less if the quarks have equal masses. In that case A = and the matrix D 
for the the unnatural parity P = (— j''^^ states falls apart into two 4x4 parts describing the charge-parity C — {—Y 
states and the C = {—)'''^^ states 

{V}, 'rj, 3(J-1)}, ^(J+l)}} and fj}, VJ, HJ-ITj, HJ+m, (All) 

respectively. If the masses are equal the set of natural parity P — (—)"■' states can be split into the C = (— states 
and the C = (— )'^^^ states 

{V}, VJ, 3(J-1)}, HJ-ITj, HJ+1):j, HJ+m and {ijj, V}}, (A12) 

respectively. The latter C = (~)'^^^ set contains no |+ +) components and therefore it has no non-relativistic analog. 
However since Eq. ( |A9|) is independent of E for this case, no bound states exist in this sector. 



2. Equal-time approximation 



As with the BSLT propagator, the parts Si{x) and S2{x) in Eqs. (2.32) and (2.33) of the ET propagator can be 
projected on the partial wave bases. One finds for the set B (A6) of unnatural parity states 

Si = mim2l + 
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dj 


—miE2 








ciAD-j+i 


ciE'2-D-j+i 


~C2ADj+2 


-C2i?2-Dj+2 \ 






-d'j 








-ciE2D^j+i 


ciBD^j+i 


C2-E2-DJ+2 


-C2BDJ+2 














C2BD_j+i 


—C2E2D-J+1 


ciBDj+2 


— C1E2DJ+2 










-TOii?2 


dj 


— C2i?2-D_j+i 


C2AD_j+i 


CIE2DJ+2 


ClADj+2 




ciADj+i 


~ciE2Dj+i 


C2BDJ+1 


C2£'2-Dj+i 


l_a2 


miE2 


-2C1C2H2 







C1E2DJ+1 


ciBDj+i 


-C2E2DJ+1 


C2ADJ+1 


miE2 


2J+l"j-l 





2C1C2H2 




-C2AD.J 


C2E2D^j 




ciE2D^j 


-2C1C2H1 





1 f)2 


miE2 
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-C2E2D^j 


-C2BD^j 


~ciE2D^j 







2ciC2i?l 


TOli?2 


2J+l"j+l / 



with Dg, ci, C2 as defined in (A3), A = mi + 7712, B = mi — 7712, and 

Hi 



H2 



92 


-2J+1 d 


-J2 + 1 




X dx 


2a;2 


92 


2J + 3d 


J2-2J 


dx^ 


X dx 


2x2 


52 


2d e(£+l 


) 


dx^ 


X dx 2x^ 





The projection of Si for the set B* (A7) of natural parity states is given by 

SI - TSiT, T = diag(l, -1, 1,-1, 1, -1, 1,-1). 
The projection S2 of ^2 for the set B can be expressed as 

S2imi,m2,Ei,E2) = CSi{m2,mi, E2, Ei)C , 
with diagonal charge-parity matrix C = diag(l, 1, —1, —1, 1, —1, 1, —1). And finally, 



(A13) 

(A14) 
(A15) 
(A16) 

(A17) 
(A18) 
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S*=TS2T or Slimi,m2,Ei,E2) = C*Sl{m2,mi,E2,Ei)C*, 



(A19) 



with charge-parity matrix C* = TC ~ diag(l, —1, —1, 1, 1, 1, 1, 1). The complete partial wave projected ET equation 
becomes the following set of eight coupled integro-differential equations 



/ 



dx X 



— Zt{x,x')Si{x') + —Z2{x,x')S2{x') 

mi 7712 



ip{x') = -V{x)xl}{x). 



(A20) 



The diagonal matrices Zi have as j-th diagonal element the projection of Zi on the orbital angular momentum tj of 
the j-th basis state, 



Zf\x, x) = 27r / d7 F£(7) Z^{y/ x"^ + x'"^ - 2xx'-f). 



(A21) 



3. Potential 



The vector structure of the qq-interaction gives 



1+f 
1-/ 









-i + f 



(A22) 



-1-/, 



on the basis (|A4|). Here / = — 3 for spin singlet states and / = 1 for spin triplet states. Choosing the Coulomb gauge 
amounts to replacing 



1 



(cr(i) •a;)(^(2) 



(2) 



(A23) 



Thus in the Coulomb gauge the spin-spin interaction is reduced. The difference with the Feynman gauge can be 
expressed as 
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-CiC2 





^2 
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(A24) 



both for the bases Mfl) and (A7). The constants Ci have been given in Eq. ([A3|). 



APPENDIX B: SHORT DISTANCE BEHAVIOR OF ^So STATES 



Here we discuss in detail the singular behavior of the wave function of the = 0~ states (e.g. pion), to which we 
refer to as states for convenience because their wave functions contain mostly ^5*0 components. At short distances, 
or equivalently at high momenta, the masses of the two quarks can be neglected, so the short distance behavior for 
the unequal quark mass case is identical to the equal quark mass case. Thus in this appendix only the equal mass case 
needs to be considered. Since the short distance behavior for the ET equation has already been discussed in Ref. Q| 
we concentrate on the BSLT approximation. Since J = the J — 1)} state is removed from the basis in ( All ) and 
the wave function is decomposed as 



(Bl) 



In this basis the BSLT equation ( |A9| ) reduces to 
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\ ■03 (a:) 



Z(o)(a;,.T') 





/ m E Di 
E m 
\ D-i m 



/-2 00\ /lOO 

4 Vv{x') + 1 I Vs{x') 
V / V 1 



(B2) 



In this formula E ~ M/2 and VV(2;) and Vs{x) denote the vector and scalar part of the qq potential. The Feynman 
gauge is assumed. We first consider the case of a fixed coupling constant a. Assuming that ipi{x) = Cix'^^ for x — > 
one can reduce in leading order in x the BSLT equation to 



— mcix''^ + Ec2X^^ + 03(73 + l)a;^-''~"'" 
Ecix'^^ — mc2x'^^ 
ciili — l)x'^^~^ — mcsx'^^ 



-2mcix''^I^°'> (71) + AEc2x'^U^°'> (72^ 
-2ci(7i-2)x^i-i/(i)(7i-l) 



(B3) 



where the /*^^^ follow from 



m .10 



dx' xx' Z^'^\x,x') x'-^-^ = x^l'^^X^), 



(B4) 



which is well-defined if a; and I7I < 1. Assuming that Eq. (B3) does not solely determine the energy eigenvalue 
E, we may exclude the case that 71 = 72. As a result we find that the exponents satisfy 72 < 71, 73 = 72 + 1 and 
= /^°''(72). Taking the explicit form of 



Z(0)(a;,a;') 



KQ(m\x - x'\) - KQ{m\x + x'\) 



we get 



/^°'(7)--r^ln(i±^) tan ^. 

Jo V \^-vJ 7 2 



Explicit expressions for general /*^^^ are given in Ref. 0]. Hence 



a 



1 2 

< 



(BSLT) 



This may be compared with the asymptotic behavior of the BSE pj| , where 



72 = 73 - 1 



a<-, (BSE) 



7i=74 = \/5-2\/4+^(l + ^ 



(BSE) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



Here 74 is the exponent of the ^Pq component which does not decouple in the BSE. Note that the dependence of 72 
on a is to first order identical for the BSLT equation as for the BSE. 

The ET uv behavior can be found in a similar way and was already discussed in Ref. Here again 72 < 71, 
73 = 72 + 1 and 



a= i72(l-72)/^°Hl-72) < -• (ET) 

Z TT 



(BIO) 
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If the Feynman gauge is replaced by the Coulomb gauge then according to Eq. (A24) we have to replace the diagonal 
matrix (—2, 4, 0) in Eq. (B2) by the diagonal matrix (—1, 3, —1) and one gets for the BSLT and ET equations 



a = ^72(l-72)/(°)(l-72)<^. (ET) (B12) 

So the short distance behavior is less singular in the Coulomb gauge than in the Feynman gauge. This can be 
understood from the suppression of the spin-spin interaction in the Coulomb gauge. The same singular behavior has 
been found by Murota Q for the Salpeter equation, which is obtained from the ET equation by removing the |H — ) 
and I — h) components. 

If the short distance attraction of the Coulomb-like potential is weakened by taking a running coupling constant 
a{x) ~ ln~^(a;) then the singular short distance behavior of the tjji is also weakened to = Ci| ln(x)|'''». However, 
to construct numerically the solutions to the quasi-potentials equations it is not necessary to know the expressions 
for the "fi. The equations can be solved by assuming the singular behavior associated with a fixed coupling constant 
which is equal to the running coupling constant at a typical small distance. 
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FIG. 1. The meson wave function at zero relative distance written as a loop integral over the quark (quasi-) propagators and 
the potential, with additional corrections. If this loop integral is divergent, then the wave function is divergent at x = 0. 

FIG. 2. A physical bound state ^So solution (a) at M = 1.23 GeV and an unphysical continuum solution (b) at 
M = 1.30 GeV obtained from the BSLT equation, with m = 0.2 GeV, k = 0.2 GeV^, a = c = e = 0. 

FIG. 3. The masses of the ground state and the first few excitations of the ^So systems are shown as a function of the fixed 
coupling constant a (calculated with m = 5 GeV, k = 0.2 GeV^, e = c = 0). The solid lines are for BSLT, broken lines for 
ET, and dotted for Schrodinger. 

FIG. 4. The masses of the ground state and the first few excitations of the ^So system are shown as a function of the constant 
term c in the qq potential (calculated with m — 0.2 GeV, k — 0.2 GeV^, e = a = 0). States with negative mass represent 
anti-mesons. The solid lines are for BSLT, broken lines for ET. 

FIG. 5. The masses of the ground state and the first few excitations of the ^5*0 systems are shown as a function of the quark 
mass (calculated with c = £ = a = 0), but with the ratio w/m^ = 0.2 fixed. Under this condition the Schrodinger levels are 
independent of m; they are indicated by the arrows on the right of the figure. The solid lines are for BSLT, broken lines for 
ET. 



FIG. 6. The masses of the ground state and the first few excitations of the ^5*0 systems in the BSLT equation are shown as 
a function of the mass mi of the first quark (calculated m2 = 0.2 GeV, k = 0.2 GeV^, c = e = a = 0). 



FIG. 7. Regge-slopes in the BSLT model. The squared meson masses as a function of the total angular momentum J of 
the unexcited mesons (calculated with m — 0.25 GeV, a = 0, c = —1.0 GeV and k = 0.33 GeV^). Solid lines show the 
masses for a pure scalar confinement (e = 0) , broken lines for mixture of vector and mainly scalar confinement in the Feynman 
gauge {e = 0.15), and dotted lines for the same mixture using the Coulomb gauge. The lines labeled n represent the vr family 
[P = (-)-'+^ and C = {-y], and similarly for the p family [P = C = (-)-'+^] and the m family [P = C = {-)■']. The 
corresponding Regge-slopes are listed in Table |. 



FIG. 8. Identical to Fig. \A but now for the ET equation. 



TABLE I. Slope f3 (in GeV'^) of the Regge-trajectories shown in Fig. ^ and Fig. ^. Listed values 
follow from a fit from J = 2 to J = 6. 





family 


e = 


e = 0.15 
Feynman 


e = 0.15 
Coulomb 


BSLT 


TT 


0.69 


1.20 


1.06 




P 


0.69 


1.13 


1.01 




ffll 


0.66 


1.12 


1.01 


ET 


TT 


0.67 


0.47 


0.64 




P 


0.69 


1.04 


0.99 




ai 


0.65 


0.92 


0.83 
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TABLE II. Parameters for the BSLT and ET model (in units of GeV). 





BSLT 


BSLT 


ET 




Feynman gauge 


Coulomb gauge 


Coulomb gauge 




0.200 


0.250 


0.250 


TTls 


0.404 


0.447 


0.390 


rUc 


1.715 


1.779 


1.719 


rub 


5.121 


5.199 


5.096 


K 


0.33 


0.33 


0.33 


c 


-0.8 


-1.0 


-1.0 


e 


0.2 


0.25 


0.2 




0.8 


0.8 


0.8 


running type 


ai, /J = 1 


ai, 11 = 1 


an 
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TABLE III. The meson mass spectrum (in GeV). The experimental data are from the particle 
data group Numbers between brackets need confirmation. Only experimental errors larger 

than 0.01 GeV are given. The column labeled N^^^^Lj displays the quantum numbers of the main 
component of the wave function. Underlined values have been fitted. The most satisfying results 
are found from the BSLT using the Feynman gauge. 

quark meson J^'^ N^'^+^Lj Observed BSLT BSLT ET 

content mass Feynman Coulomb Coulomb 



TT 


X 

+ 


z — n^-, 

1 'So 


0.135 


0.439 


0.521 


0.134 


tt' 


0"+ 


2 ^Sn 


1 30+ 10 


1.441 


1.424 


1.361 


ir" 


0-+ 


3 'So 


(1.77± 0.03) 


2.246 


2.193 


1.978 


P 


1 


1 ^Si 


0.768 


0.798 


0.796 


0.838 


P' 


1 


1 ^Di 


1.47± 0.03 


1.454 


1.470 


1.613 


P" 


1 


2 ^Si 


1.70± 0.02 


1.653 


1.594 


1.649 


o'" 
P 


1~~ 


2^Di 


(2 14+ 03"! 


2.185 


2.180 


2.230 


P"" 


1 — 


3 -'Si 




2.367 


2.308 


2.248 


bi 


1+ 


1 Pi 


1.23 


1.091 


1.136 


1.196 


ao 


0++ 


1 ^Po 


0.983 


0.993 


1.001 


1.186 


n.-\ 


1++ 


1 


1 26± 03 


1.126 


1.142 


1.249 


02 


2++ 


1 ^P2 


1.318 


1.297 


1.277 


1.311 


772 


2 + 


1 '-D2 


1.67± 0.02 


1.524 


1.552 


1.569 


63 


3+" 


1 




1.886 


1.896 


1.820 


71'4 


4"+ 


1 




2.205 


2.198 


2.011 




5+- 


1 ^Hr 

_L i J b 




2.494 


2.471 


2.164 




6-+ 


1^/6 




2.761 


2.722 


2.236 




3 


1 ^Ds 


1.69 


1.689 


1.654 


1.660 




4++ 


1 ^ Pi 


(2 04+ 03"! 


2.021 


1.973 


1.946 




5-- 


1 ^Gs 


(2.35) 


2.314 


2.253 


2.191 


oe 


6++ 




(2.45+ 0.13) 


2.583 


2.507 


2.407 


P2 


2— 


1 ^D2 




1.547 


1.562 


1.624 


as 


3++ 


I'Ps 


(2.08+ 0.04) 


1.897 


1.903 


1.904 


PA 


4-- 


1 ^G4 




2.202 


2.199 


2.133 


as 


5++ 


1 ^^5 




2.479 


2.465 


2.330 


P6 


6-~ 


l'/6 




2.732 


2.709 


2.503 


K 


0" 


1 'So 


0.498 


0.593 


0.660 


0.350 


K' 


0" 


2 'So 


(1.46) 


1.560 


1.530 


1.457 


K" 


0" 


3 'So 


(1.83) 


2.326 


2.268 


2.069 


K* 


1" 


1 ^Si 


0.896 


0.896 


0.896 


0.910 


K*' 


1" 


1 ^Di 


1.41 


1.584 


1.599 


1.692 


K*" 


1" 


2 ^Si 


1.71+0.02 


1.727 


1.680 


1.700 


Ki 


1+ 


1 'Pi 


1.27 


1.232 


1.257 


1.282 


K*o 


0+ 


I'Po 


1.43 


1.113 


1.129 


1.262 


Ki 


1+ 


1 ^Pl 


1.40 


1.261 


1.268 


1.333 


Kl 


2+ 


1 ^P2 


1.425 


1.389 


1.373 


1.376 


K2 


2" 


1 'D2 


1.77 


1.655 


1.669 


1.672 


Kz 


3+ 


1 'P3 


(2.32+0.02) 


2.005 


2.008 


1.943 


Ki 


4- 


1 'G4 


(2.49+0.02) 


2.331 


2.304 


2.153 


K5 


5+ 


1 'H3 




2.590 


2.571 


2.326 


Kl 


3" 


1 ^Ds 


1.77 


1.776 


1.749 


1.726 


Kl 


4+ 


1 ^Pi 


2.05 


2.107 


2.069 


2.018 


Kt 


5- 


I'Gs 


(2.38+0.03) 


2.399 


2.350 


2.271 


D 


0" 


1 'So 


1.864 


1.868 


1.912 


1.855 


D" 


1" 


1 ^Si 


2.007 


2.015 


2.032 


2.058 
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Di 




1 1 

1 Pi 


u 


0+ 


1 ^Po 




1+ 


1 


Di 


2+ 


1 


Ds 


0~ 


1 1 rr 

1 ^5o 




1" 


1 ^Si 


L>s\ 


1 + 
i 


1 Pi 


Dto 


0+ 


1 ^Pl 


D'si 


1+ 


1 ^Pl 




2+ 


1 ^P2 


Vc 


U 


1 1 c 

1 00 


Vc 


n — h 


1 

2 Ao 


J/iJ 


1 — 




V'' 


1 — 




tp 


1 


-1 3 n 

1 l?! 


1 fff 


1 


3 


1 iv 


1 


2 


V 


1 — 




hci 




1 1 
1 Pi 




0++ 


1 ^Po 


Xcl 


1++ 


1 ^Pi 


Xc2 


2++ 


1 3P2 


B 


0~ 


1 1 rr 
1 So 


B* 


1" 




Bi 




1 1 
1 Pi 


Bn 


0+ 


1 ^Po 


B[ 


1+ 


1 3Pl 




2+ 


1 ^P2 


Bs 


0~ 


lie 
1 00 


B*s 


1" 


1^51 


Bsi 


1 + 


1 1 D 
1 Pi 


B^n 


0+ 


1 ^Po 


Bsi 


1+ 


1 ^Pl 


B32 


2+ 


1 3P2 


Be 


0~ 


lie 

1 60 


B* 


1" 






Del 


1 + 
i 


1 Pi 


0* 
JDcQ 


0+ 


1 ^Pi 




1+ 


1 ^Pi 


B*2 


2+ 


1 ^P2 


Vb 


^ 


lie 

1 So 


Vb 


^ 


2 ^So 


T 


1 — 


l^Si 


^/ 


1 — 


2 ='S'i 


T 


1 


1 ^Di 


1 


1 

1 


Q 3q 
Dl 


^iv 


1 — 


2^P)i 


'^v 


1 — 


4351 


hbl 


1+- 


l^Pi 


Xbo 


0++ 


l^Po 



2.42 


2.388 


2.404 


2.423 




2.321 


2.327 


2.632 




2.415 


2.420 


2.494 


2.459 


2.458 


2.461 


2.474 


1.969 


1.952 


1.983 


1.939 


(2.110) 


2.104 


2.110 


2.115 


2.537 


2.500 


2.509 


2.491 




2.427 


2.419 


2.510 




2.516 


2.502 


2.550 


(2.564) 


2.569 


2.559 


2.547 


2.980 


2.969 


2.988 


2.883 


(3.59) 


3.742 


3.713 


3.680 


3.097 


3.096 


3.098 


3.097 


3.686 


3.810 


3.779 


3.789 


3.770 


3.873 


3.854 


3.872 


4.04 


4.370 


4.325 


4.308 


4.16± 0.02 


4.409 


4.372 


4.363 


4.42 


4.860 


4.805 


4.744 


3.526*^ 


3.517 


3.513 


3.486 


3.415 


3.461 


3.442 


3.486 


3.511 


3.526 


3.510 


3.521 


3.556 


3.572 


3.559 


3.553 


5.279 


5.302 


5.331 


5.349 


5.324 


5.360 


5.383 


5.391 




5.741 


5.756 


5.741 




5.714 


5.717 


5.800 




5.760 


5.766 


5.817 




5.770 


5.781 


5.760 


(5.38±0.03) 


5.371 


5.383 


5.390 


4.3+n 03^ 


5.434 


5.443 


5.441 




5.839 


5.841 


5.805 




5.802 


5.789 


5.848 




5.846 


5.838 


5.865 




5.869 


5.866 


5.827 




6.260 


6.260 


6.228 




6.331 


6.329 


6.336 




6.760 


6.754 


6.719 




6.724 


6.702 


6.724 




6.767 


6.745 


6.742 




6.794 


6.781 


6.758 




9.401 


9.402 


9.350 




10.067 


10.047 


9.980 


9.460 


9.460 


9.459 


9.460 


10.023 


10.099 


10.081 


10.039 


10.355 


10.206 


10.187 


10.138 


10.58 


10.556 


10.532 


10.467 


10.87 


10.629 


10.603 


10.533 


11.02 


10.943 


10.911 


10.827 




9.881 


9.879 


9.823 


9.860 


9.862 


9.843 


9.831 
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Xbl 


1++ 


1 


9.892 


9.890 


9.876 


9.845 


Xb2 


2++ 


1 ^P2 


9.913 


9.911 


9.901 


9.863 


Kl 


1+- 


2 ^Pi 




10.383 


10.363 


10.288 


XbO 


0++ 


2'Po 


10.232 


10.363 


10.335 


10.289 


X'bl 


1++ 


2 ^Pi 


10.255 


10.384 


10.360 


10.301 


X'b2 


2++ 


2^P2 


10.268 


10.400 


10.379 


10.316 



" Ref. 1^. 



TABLE IV. 16-component vector qq-basis states and cor- 
responding 4x4 matrix qq-basis states. The angular de- 
pendence of the wave functions of the qq-triplet states must 
determine to which qq-triplet state they correspond. For a 
correct normalization all matrix states should be multiplied 
by 1/2; furthermore, states proportional to P or p need an 
extra factor 1/M and l/|p|, respectively. The correspondence 
is only valid in the cm system and po = is assumed. 



qq-state qq-state 

^Jj 75 

V5 -1 

V}, Hj-iYj, HJ+m - 

HJ-ITj, ^(J+1)} ia^.P^'p- 

^J^j, HJ-ITj, 3(J+1)} -75ia^.PV 

'r,, '{J-m, ^(J+1)} 75 
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